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We present the results of a LIGO search for short-duration gravitational waves (GWs) associated 
with Soft Gamma Repeater (SGR) bursts. This is the first search sensitive to neutron star /-modes, 
usually considered the most efficient GW emitting modes. We find no evidence of GWs associated 
with any SGR burst in a sample consisting of the 27 Dec. 2004 giant fiare from SGR 1806—20 and 
190 lesser events from SGR 1806—20 and SGR 1900-1-14 which occurred during the first year of 
LIGO's fifth science run. GW strain upper limits and model-dependent GW emission energy upper 
limits are estimated for individual bursts using a variety of simulated waveforms. The unprecedented 
sensitivity of the detectors allows us to set the most stringent limits on transient GW amplitudes 
published to date. We find upper limit estimates on the model-dependent isotropic GW emission 
energies (at a nominal distance of 10 kpc) between 3 x 10*^ and 9 x 10^^ erg depending on waveform 
type, detector antenna factors and noise characteristics at the time of the burst. These upper limits 
are within the theoretically predicted range of some SGR models. 



Soft Gamma Repeaters (SGRs) sporadically emit brief 
(rs 0.1s) intense bursts of soft gamma-rays with peak 
luminosities commonly up to 10^^ erg/s Q. Less 
common intermediate bursts with greater peak luminosi- 
ties can last for seconds. Rare "giant flare" events, some 
1000 times brighter than common bursts have initial 
bright, short (« 0.2 s) pulses followed by tails lasting min- 
utes and are among the most electromagnetically lumi- 
nous events in the Universe Q . Since the discovery of 
SGRs in 1979 three of the four confirmed SGRs have 
produced a giant flare each [1, [^, @] . 

SGRs are promising sources of gravitational waves 
(GWs). According to the "magnetar" model SGRs are 
neutron stars with exceptionally strong magnetic fields 
~ 10^^ G Q- SGR bursts may result from the interaction 
of the star's magnetic field with its solid crust, leading to 
crustal deformations and occasional catastrophic crack- 
ing 0] with subsequent excitation of the star's nonra- 
dial modes l3, nil, nj and the emission of GWs [nl, [H, 
Excitation of nonradial modes could also occur if 
SGRs are instead solid quark stars isl. 14. isj. 

We present a search for short-duration GW signals 
(< 0.3 s) associated with SGR bursts using data collected 
by the Laser Interferometer Gravitational Wave Obser- 
vatory (LIGO) . LIGO consists of two co-located 
GW detectors at Hanford, WA with baselines of 4 km 
and 2 km and one 4 km detector at Livingston, LA. GW 
data from one or two of these detectors are used. When 
three detectors are operating, data from the most sensi- 
tive pair are chosen. 

The SGR burst sample was provided by gamma-ray 
satellites of the interplanetary network [l3l , and includes 
the 27 Dec. 2004 giant flare from SGR 1806-20 and 
214 conflrmed bursts (152 from SGR 1806-20 and 62 
from SGR 1900-1-14, one of which was a multi-episodic 
"storm" [3]) occurring during the first year of LIGO's 
fifth science run (S5) from 14 Nov. 2005 to 14 Nov. 2006. 



Of the 214 bursts, 117 occurred with three LIGO detec- 
tors operating, 53 with two detectors operating, 20 with 
a single detector operating, and 24 with no detector op- 
erating. Including the giant flare, analysis was possible 
for a total of 191 listed SGR events. 

To analyze a given SGR burst we divide the GW data 
into an on-source time region, in which GWs associated 
with the burst could be expected, and a background time 
region. In the background region we do not expect a GW 
associated with the SGR burst, but the noise is statisti- 
cally similar to the on-source region. For isolated bursts 
the on-source region consists of 4 s of data centered on the 
SGR burst. GW emission is expected to occur almost 
simultaneously with the electromagnetic burst [l^ ; the 
4 s on-source duration accounts for uncertainties in the 
geocentric electromagnetic peak time. There are three 
special cases: 1) for two SGR 1900-1-14 bursts which oc- 
curred within 4 s of each other a combined 7 s on-source 
region was chosen; 2) for the SGR 1900+14 storm a 40 s 
on-source region was used; 3) for the GRB 060806 event 
from SGR 1806—20, two 4 s on-source regions were used, 
centered on the two distinct bright bursts comprising the 
event. Background regions consist of 1000 s of good data 
on either side of on-source regions. On-source and back- 
ground segments are analyzed identically, including data 
quality cuts, resulting in lists of "analysis events." Anal- 
ysis events from the background regions are used to es- 
timate the significance of the on-source analysis events; 
significant events, if any, are subject to environmental 
vetoes and consistency checks. 

The analysis is performed by the Flare pipeline 3, 2^, 



2l| and is based on the excess power detection statistic 
of j2^ . Search parameters such as frequency bands and 
time windows are chosen to optimally detect the target 
signals. This is achieved by comparing detection efficien- 
cies for simulated target signals injected into the back- 
ground data and searched for with different search pa- 
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rameters [20| . The search targets neutron star fundamen- 
tal mode rine; downs (RDs) predicted in [l^, [O, [H, El, [H 
as well as unmodeled short-duration GW signals. Model 
predictions from pjj for ten realistic neutron star equa- 
tions of state give /-mode RD frequencies in the range 
1.5-3 kHz and damping times in the range 100-400 ms. 
We use a search band 1-3 kHz for RD searches (to in- 
clude stiffer equations of state), and find a 250ms time 
window to be optimal. The search for unmodeled signals 
uses time windows set by prompt SGR burst timescales 
(5-200 ms) and frequency bands set by the detector's sen- 
sitivity; a 125 ms time window effectively covers this du- 
ration range, and we search in two bands: 100-200 Hz 
(probing the region in which the detectors are most sen- 
sitive) and 100-1000 Hz (for full spectral coverage below 
the ringdown search band). 

In the absence of a detection, for each SGR burst 
we estimate loudest event upper limits [2^ on the GW 
strain incident on the detector, h^-ss- Following 



/?2 

"tss-I- 



where e.g. 



"rss+ 



h^dt and 



{t) are the two GW polarizations. The relationship 
between the GW polarizations and the detector response 
h(t) to an impinging GW from an altitude and azimuth 
{6, (p) and with polarization angle V' is: 

h{t) ^ F+{0,<f>,^)h+{t) + F^{0,cp,^)h^{t) (1) 

where F^{9, 0, ip) and (9, 0, ip) are the antenna func- 
tions for the source at {9, (j)) j28|. The upper limit is com- 
puted in a frequentist framework following the commonly 
used procedure of injecting simulated signals in the data 
and recover ing them using the search pipeline (see for ex- 
ample (29l . l30l|). The upper limits are derived for RD 
signals and for unmodeled bursts. Correspondingly, RD 
and band-limited white noise burst (WNB) waveforms 
are injected with parameters chosen to probe the respec- 
tive target signal space. For WNBs independent polar- 
ization components are generated with /iiss+ = ^rssx- 
For RD signals linearly and circularly polarized waves 
are considered. 

The GW strain h-^ss upper limits can be recast as up- 
per limits on the emitted GW energy, -Eqw- Assuming 
isotropic emission, the GW energy associated with /i+ [t) 
and /ix {t) is [3ll |: 



GW 



IQttG 



({h+f + {h^f^ dt. (2) 



We use this equation with a nominal source distance of 
R = lOkpc (source locations and distances are discussed 
in (3^. [33j ) to compute the energies associated with the 
/irss upper limits for different signals. 

Results — We find no evidence for gravitational waves 
associated with any of the SGR burst events in the sam- 
ple. The significance of on-source analysis events is in- 
ferred by assigning rates at which background analysis 
events of equal or greater loudness occur. Wc find the 



most significant on-source analysis event occurs at a rate 
of 1.35 X 10"'^ Hz (1 per 741s), which is consistent with 
the expectation for the 803 s of on-source data in the 
sample. We estimate 90% confidence strain and energy 
upper limits, h^ss° ^^'i -^gW' using the loudest on-source 
analysis event for each SGR burst. Upper limits depend 
on detector sensitivity and antenna factors at the time of 
the burst, the loudest on-source analysis event, and the 
simulation waveform type used. 

TableU hsts upper limits for the SGR 1806-20 gi- 
ant flare and for the brightest peak of the GRB 060806 
event (33| (complete results are given in 35|). Results 
from these two events are highlighted because they yield 
the smallest values of 7 = E'^^/E-em^ a measure of 
the extent to which an energy upper limit probes the 
GW emission efficiency. At the time of the giant flare 
the LIGO Hanford 4 km detector was operating during a 
commissioning period (LIGO Astrowatch) and had noise 
amplitude higher than that of S5 by a factor of 3; the 
rms antenna factor, which is an indicator of the aver- 
age sensitivity to a given source in the sky, for such 
event was (Fj -l-F^)^/^ = 0.3. The isotropic electromag- 
netic energy {E-em) for the event, assuming a distance 
of lOkpc, was 1.6 x lO^'^ergjil. At the time of GRB 
060806 both the 4 km and 2 km Hanford detectors were 
observing, with rms antenna factor for that event of 0.5. 
ii^EM for the brightest peak of GRB 060806 was at least 
2.9 X 10*2 ergjsj. 

We estimate upper limits on GW strain and isotropic 
GW energy emitted using RDs with t = 200 ms at var- 
ious frequencies, and WNBs lasting 11 and 100 ms and 
with 100-200 and 100-1000 Hz bands. We observe no 
more than 15% degradation in strain upper limits us- 
ing RDs with r in the range 100-300 ms, and no more 
than 20% degradation using WNBs with durations in 
the range 5-200 ms, as compared to the upper limits ob- 
tained for the nominal RDs and WNBs used for tuning 
the search. Superscripts in Table|T] give a systematic er- 
ror and uncertainties at 90% confidence. The first and 
second superscripts account for systematic error and sta- 
tistical uncertainty in amplitude and phase of the detec- 
tor calibrations, estimated via Monte Carlo simulations, 
respectively. The third is a statistical uncertainty aris- 
ing from using a finite number of injected simulations, 
estimated with the bootstrap method using 200 ensem- 
bles [3^ . The systematic error and the quadrature sum of 
the statistical uncertainties are added to the final upper 
limit estimates. 

Figure[T] shows E'^^ limits for the entire SGR burst 
sample. The lowest upper limit in the sample, E"^^ = 
2.9 X lO""^ erg, is obtained for an SGR 1806-20 burst on 
21 Jul. 2006, with a geocentric crossing time of 17:10:56.6 
UTC. The lowest upper limit from the RD search is 
^90% = 2.4 X lO'^^erg for an SGR 1806-20 burst on 
24 Aug. 2006 14:55:26 UTC. 

Discussion — Two searches for GWs associated with 
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FIG. 1: E^-^ upper limits for the entire SGR burst sample 
for various circularly /linearly polarized RD (RDC/RDL) and 
white noise burst (WNB) signals. The limits shown in TableU] 
for the giant flare and GRB 060806, are indicated in the figure 
by circles and diamonds, respectively. 
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SGR events have been published previously; neither 
claimed detection. The AURIGA collaboration searched 
for GW bursts associated with the SGR 1806—20 giant 
flare in the band 850-950 Hz with damping time 100 ms, 
setting upper limits on the GW energy of ~ 10^^ erg [sj- 
The LIGO collaboration also published on the same gi- 
ant flare, targeting times and frequencies of the quasi- 
periodic oscillations in the flare's x-ray tail as well as 
other frequencies in the detector's band, and setting up- 
per limits on GW energy as low as 8 x 10''^ erg for quasi- 
periodic signals lasting tens of seconds [sl] . 

In addition to the 2004 giant flare, the search described 
here covers 190 lesser events which occurred during the 
LIGO S5 data run. Furthermore this search extends to 
the entire high sensitivity band of the detectors, which 
makes it the first search sensitive to neutron star /- 
modes, usually considered the most efficient GW emitting 
modes [l^l . Our upper limits on i?GW overlap the range 
o{Eem lO'^^'-lO'^^ erg seen in SGR giant flares [l,!!]. Most 
of the WNB limits, and some of the RD limits, are below 
the 10*^^ ergmaximum E'gw predicted in some theoret- 
ical models |12| . Our best upper limits on t are within 



the theoretically predicted range implied in | 

The Advanced LIGO detectors promise an improve- 
ment in /irss by more than a factor of 10 over S5, corre- 
sponding to an improvement in energy sensitivity (and 
therefore 7) by more than a factor of 100. Thus within 
the next few years we expect to obtain GW energy upper 
limits for the /-mode search that fall in the Eem range 



of giant flares, and for the unmodeled search that fall in 
the -Eem range of intermediate bursts. 
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TABLE I: GW strain and energy upper limit estimates at 90% confidence {h^ss and -Eqw) SGR 1806—20 giant flare and 

the S5 SGR burst with the smallest limits on the ratio 7 = -Eq^Z-Eem for various circularly /linearly polarized RD (RDC/RDL) 
and white noise burst (WNB) waveforms. Uncertainties (given in superscripts for strain upper limits and explained in the text) 
are folded into the final upper limit estimates. 
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